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The mutual influence between a Kondo impurity and a d-wave superconductor is numerically studied 
within the slave-boson mean-field approach. The Bogoliubov-de Gennes type equations are derived 
and solved using exact diagonalization. We show that a critical coupling strength, above which the 
Kondo efi^ect takes place, exists regardless of whether the band particle-hole symmetry is present or 
not. In the Kondo regime the double resonant quasiparticle peaks are found in the local density of 
states (LDOS) both directly at the impurity and around its neighbors, which is in sharp contrast to 
the case of nonmagnetic unitary impurities, where the LDOS vanishes at the impurity site. 
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The identification of c?-wave pairing state in high-Tc 
superconductors has renewed much interest in low- 
energy quasiparticle states in unconventional supercon- 
ductors. Different from conventional s-wave supercon- 
ductors having a constant energy gap, the d-wave super- 
[ conductors have an anisotropic energy gap which changes 
sign along the nodal directions on the essentially cylin- 
drical Fermi surface. The anisotropy of the energy gap 
makes the d-wave superconductors very sensitive to im- 
purity or other potential scattering, in contrast to the 
s-wave superconductors which, as dictated by the An- 
derson theorem |^], are almost unaffected by the pres- 
ence of nonmagnetic impurities. As an important sig- 
nature of the sign change of d-wave order parameter 
(OP), Balatsky, Salkola and co-workers based on 

the T-matrix approximation, theoretically predicted that 
a strong atomic-like nonmagnetic impurity can induce 
a virtual bound state. Later studies showed that 
the resonant state would show up in the local density 
of states (LDOS), which could be measured by the local 
differential tunneling conductance. Recent scanning tun- 
neling spectroscopy measured on high-T^ superconduc- 
tor Bi2Sr2CaCu208+5 (BSCCO) with individual atomic- 
scale defects or strong impurities has provided 
strong evidence for the existence of low-energy quasi- 
bound states. Especially, the experiment by Pan et al. ||] 
on the observation of quasiparticle resonant states has a 
high spatial and energy resolution and the identity of 
impurity atoms is known in a controlled way. Very inter- 
estingly, the imaging of the effect of a single zinc (Zn) im- 
purity atom showed the strongest intensity of the LDOS 
at resonance directly on the Zn atom, which is opposite 
to the earlier theories . In our previous work [|o| , we 
explained the novel pattern of the tunneling conductance 
in terms of the blocking effect of the BiO and SrO lay- 
ers which exist between the tunneling tip and the Cu02 
layer being probed. This explanation may not exclude 
other possibilities and needs further experimental justi- 
fication |l^. Although almost all existing theories are 
treating the Zn atom as a static potential scatterer, there 
exists a few experiments n% indicating that Zn impurity 



may possess a magnetic moment. It is now an appropri- 
ate time to consider the effects of dynamic impurities on 
the quasiparticle properties of d-wave superconductors. 
The study of a dynamic impurity in a superconductor be- 
gan with the works of Yu jl^ and Shiba on magnetic 
impurities in s-wave systems. The solving of this problem 
is quite challenging mainly due to the difficulty in treat- 
ing the dynamical correlations of the coupled impurity- 
conduction electron system with pair correlation. The 
problem of Kondo impurities in a d-wave system has 
been studied by focusing on either the moment forma- 
tion [^5],^ or the suppression of superconducting tran- 
sition temperature by them JlTf . Quite recently, within 
a Kondo spin model, Tsuchiura et al. | [T^ ] have studied 
the LDOS around a magnetic impurity by treating the 
impurity spin as an effective local magnetic field and the 
screening of the impurity spin by the quasiparticles was 
not considered there. However, other new theories 20 
argued that the moment may be screened deep in a d- 
wave superconducting state. Experimentally, the spe- 
cific heat measurement of Zn-doped YBa2Cu3 6.95 pl| ] 
provided a possible evidence of Kondo screening in the 
superconducting state. 

In this paper we study the Kondo effects of an Ander- 
son impurity coupled with the conduction electrons of a 
d-wave superconductor. Within a slave-boson mean-field 
approach, we are able to exactly solve the composite sys- 
tem with the self-consistency for the d-wave OP included. 
It is shown that a critical coupling strength, above which 
the Kondo effect emerges, exists regardless of the band 
particle- hole (PH) symmetry being present or not. We 
also show in the Kondo regime that a double resonant- 
peak structure shows up in the LDOS both at the im- 
purity and at its nearest neighboring sites. This result 
is different from those for the static impurities that the 
LDOS at the impurity site exhibits only a single peak be- 
low (or above) the Fermi energy for a weak repulsive (or 
attractive) scattering center while vanishes in the unitary 
scattering limit. 

To describe the impurity and its interaction with the 
band electrons of the d-wave superconductor defined on 
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a two-dimensional lattice, we start with an Anderson-like 
Hamiltonian, which consists of three parts: H = Hs + 
Himp + Hs-imp- Here Hs = - Y^ii^iH + /^'^ij)cLcj-T + 
Eij(Aij4c]^ + ^yCjiCiT) + Ei^j is the Hamilto- 

nian for the host d-wave superconductor, where C;^ (cio-) 
are the creation (annihilation) operators of a conduction 
electron with spin a at the ith site, is the hopping inte- 
gral between two nearest neighbor sites, is the chemical 
potential, gij represents the strength of d-wave pairing in- 
teraction, and Ajj is the bond pairing amplitude which is 
in turn determined self-consistently as Ay = <?ij (ci|Cj|); 
Himp = Ecr ^dd],dcr -\- Urid-^ndi is the Hamiltonian for the 
Anderson impurity, where [da] is the creation (anni- 
hilation) operator of the impurity electron with spin cr, 
ndcr = dl,da, Cd is the bare impurity energy level with re- 
spect to the chemical potential of band electrons, and U is 
the on-site repulsive interaction of the impurity electrons; 
and Hs-imp = Vj^ai^l^^^s <T^'^+'^3f'^i=o+^,cr) represents 
the hybridization of the impurity with the conduction 
electrons, where without loss of generality, we have cho- 
sen the impurity to be located at the origin i = 0, 5 
are the unit vectors of the square lattice, and V is the 
coupling strength. For our purpose of studying the low 
temperature physics, we assume that the on-site inter- 
action U is infinite so that the double occupancy on the 
impurity is forbidden, which allows us to apply the slave- 
boson mean-field theory |2^. This approach has been 
applied successfully to the low-temperature properties of 
a Kondo impurity in the presence of normal state con- 
duction electrons and to transport properties of a 
quantum dot connected to a normal and to a (conven- 
tional) superconducting lead p^ , p5t , or embedded into 
a (conventional) Josephson junction pGf . By writing 
the impurity electron operator as d„ ~ b^fa, where fa 
and b are the spinon and holon operators and subject 
to the single occupancy constraint J2a fafcr + b''b = 1, 
in the mean-filed approximation, we are able to obtain 
Himp = J2a('^d + Ao)/]/ff + Ao(5o - 1), and Hs-imp 
"^'^0 E<T(Ci^o+^,<T'^'" + fyi=o+6.a)- Here the holon oper- 
ator b has been replaced by a c-number 60, Aq is the La- 
grange multiplier introduced to enforce the single occu- 
pancy constraint. Then the composite Hamiltonian can 
be diagonahzed by performing the canonical transforma- 
tion for both the conduction and the impurity electrons, 
which leads to the Bogoliubov-de Gennes equations: 

i:(5^)(|)^..(|). a. 

where {uV-, yi^y^''^^''^spose jg quasiparticle wavefunction 
with eigenvalue En, 

Hij = -^ij ~ M'^i#0,j + i'^d + Ao)^iO(5ij , (2) 

with tij ~ Vbo for (ij) = ((50) or (OS) and t otherwise, 
Ay and b^ are subject to the self-consistent condition, 



Aij = ^ Y.(ufvf* + ufvn tanh{En/2kBT) , (3) 

n 

and 

bl = 1 ^2j2{\u-\'f{E,,) + I^Sni - f{En)]} , (4) 

n 

with f{E) = [l + exp{E/kBT)]^'^ the Fermi distribution 
function. Here we have assumed the conduction and the 
impurity electrons are of different nature and take the 
pairing interaction gij = between them, and otherwise 
(7ij = g to be a constant. Notice that the eigenvalues and 
the corresponding eigenfunction amplitudes in Eq. (Q) 
is implicitly Ao-dependent. The single occupancy con- 
straint, as given by Eq. (^, is not sufficient to determine 
the values of 60 and Aq. For this purpose, we need to 
minimize the free energy, which is derived for the first 
time in terms of the quasiparticle wavefunction as: 

F - 2fcBr^ ln[l - f{En)] - J2J2[Aiiuf*vr + c.c] 

n n iJ 

+ Yl Hit-^ji + + Ao)5io5ij - /i'5i5^oj]t'>r* + C.C.} 
" iJ 

+ E^ + ^o(^o-l). (5) 

This minimization will give two equations for Aq and 
60, one of which is just the single occupancy constraint 
Eq. (jij). When applied to other physical problems, this 
expression for the free energy is very useful in determin- 
ing the stability of competing pairing states. 

Equation (^ always has a trivial solution with Aq = 
— and b^ — 0, which describes the impurity decoupled 
with the conduction electrons. This trivial critical point 
can be shown not necessarily to have the lowest free en- 
ergy. We solve Eq. (|l]) by taking Aq as a fixed parameter, 
and calculate the free energy with the obtained value of 
60 and Ay . By varying the value of Aq and repeating the 
same procedure, we obtain the free energy as a function 
of Ao . Comparing the minimum free energy with that of 
Ao = —Ed and 60 = 0, the true solution for Ao and 60 is 
then determined. The numerical calculation is performed 
on a square lattice with size Nl = Nx x Ny = 35 x 35 
and averaged over 36 wavevectors in the supercell Bril- 
louin zone. Throughout the work, we are interested in 
the zero temperature limit, and take the pairing inter- 
action g = 2t, the bare impurity level = —2t. With 
respect to the real materials, this value of g leads to an 
exaggerated amplitude of the energy gap. Our choice of 
the parameters is motivated by the desire to have a short 
coherence length. The calculations with smaller values 
of energy gap have shown similar results and no qualita- 
tive changes appear. In Fig. |l| we plot 5o and the free 
energy per site as a function of Aq for various values of 
the coupling strength V. Also shown is the value of the 
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FIG. 1. The parameter bo (a) and the free energy per site 
as a function of the Lagrange multiplier Ao for the coupling 
strength V — l.Ot (solid line), 0.8t (dashed line), and 0.6t 
(dotted line). Also shown in (b) the value of F/NLt with 
Ao — 2t and bo = Q (filled circle). The chemical potential 
/i = — 0.5t, the bare impurity level ed = — 2i, the d-wave 
pairing interaction g — 2t. 

free energy Fq obtained with Aq = — and bo = 0, which 
is independent of V. We find a critical value Vc ~ OM 
for both /i = —0.5 and 0, below which the minimum 
free energy is larger than Fq. In this regime, the im- 
purity is decoupled with conduction band, and the elec- 
tron on the impurity site is described by the spinon only. 
When V > Vc, this free local spin state is unstable against 
the couphng to the conduction electrons. We find: For 
the couphng strength V = l.Ot, 0.8t, and 0.6t, respec- 
tively, (Ao,5o) = (2.3i,0.56), (2. It, 0.50), and (2t,0) for 
fi = — 0.5t, which breaks the band PH symmetry; while 
(Ao, 5o) = (2.4t, 0.59), (2. It, 0.42), and (2.0t, 0) for ^ = 0, 
which conserves the band PH symmetry. Our numerical 
result is consistent with the conclusion by Withoff and 
Fradkin based on the Kondo exchange model that if 
the DOS of the conduction band vanishes linearly at the 
Fermi energy, the Kondo effect is suppressed for small 
values of the coupling constant. However, for the special 
case where the band PH symmetry is present, the criti- 
cal coupling strength is still found to be finite. At a first 
glance, our result seems to differ from the earlier pre- 
diction based on the Kondo spin model by Ingersent ||2|] 
that for the PH symmetry case the Kondo effect is absent 
for all values of coupling. The discrepancy with Ref. |2|] 
may come from the facts that: (i) in the present model, 
the local PH symmetry is still broken due to the asym- 
metric Anderson impurity and/or the single-occupancy 
constraint; (ii) the existence of the Kondo effect is sen- 
sitive to the local character of conduction electrons near 
the impurity; and the magnetic property of the impu- 
rity and the quasiparticle properties of superconductors 



interplay with each other, that is, the presence of the im- 
purity enhances the LDOS of the c?-wave superconductor 
(as shown below), which in turn screens the impurity 
spin. 
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FIG. 2. The local density of states directly at the impurity 
site (a) and its nearest neighbor (b) for the coupling strength 
V — It (solid line), 0.8t (dashed line), and 0.6t (dotted line). 
Also shown in (b) the bulk density of states (dash-dotted line). 
The chemical potential fi = — 0.5t. 

Once the parameter values of Aq and bo, and the OP 
are obtained, the local density of states is then calculated 
according to 



n\2 



SiEn-E) + \v]'\'S{E„ + E)] 



(6) 



which is proportional to the local differential tunneling 
conductance as measured by the STM. Figure 2 shows 
the LDOS as a function of energy at the impurity site 
and its nearest neighbor with /i — — 0.5t. Also shown in 
Fig. |2|(b) is the bulk density of states in the absence of the 
impurity. In the bulk DOS, the superconducting coherent 
peaks are exhibited at the gap edge ±Ag(= 0.75t). The 
peaks outside the gap comes from the Van Hove singu- 
larity. Due to the finite magnitude of the PH symmetry 
breaking by /i 7^ 0, the spectrum both inside and outside 
the gap is asymmetric with respect to i? = 0. More inter- 
estingly, within the energy gap, the resonant peaks show 
up in the LDOS both at the impurity site and its nearest 
neighbors, in sharp contrast to the unitary nonmagnetic 
impurity case where the LDOS vanishes at the impurity 
site The appearance of the double-peak structure 

in the LDOS directly at the impurity site also makes the 
result different from the weak nonmagnetic impurity case 
where only one peak appears below or above the Fermi 
energy {E = 0) when the scattering is repulsive or at- 
tractive 1^,0. The difference is mainly due to the fact 
that, in addition to a weak nonmagnetic-like impurity 
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scattering in the Kondo resonance phase regime, there 
is also a disruption of the effective hopping as given by 
Vbo (from the point of view of the conduction electrons) , 
immediately connecting the impurity site. The decrease 
of V enhances the quasiparticle scattering near the im- 
purity by the suppression of Vbo on one hand and shifts 
the impurity level to Cd = + -^o to the Fermi surface on 
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FIG. 3. Same as Fig. | but with = 0. 

the other hand. As a consequence of this decrease: (i) At 
the impurity site, the Kondo resonance peaks are sharp- 
ened with their separation decreased, and finally the dou- 
ble peaks merge into a single peak as V approaches Vc 
at which the decoupling of the impurity from the con- 
duction band takes place; (ii) At the nearest neighbor to 
the impurity site, the separation between double peaks 
is also decreased but the double-peak structure still sur- 
vives when V — Vc, which is intrinsic to a d-wave super- 
conductor with a short coherence when the PH symmetry 
is broken to some extent When the band PH sym- 
metry is conserved (/i = 0), as shown in Fig. ^, we see 
that most features of the local electronic structure are 
similar to the case with ^ ^ except that the spectrum 
is now symmetric outside the gap (but still asymmetric 
inside the gap due to the breaking of the local PH sym- 
metry) and that at the nearest neighbor to the impurity 
only exhibits a single zero-energy peak when V = Vc- 
We have also studied the case with different values of 
impurity level and found similar results. 

Two remarks are in order: (1) Although the slave- 
boson mean field theory captures the essential physics in 
the Kondo resonance regime, where a spin-singlet state is 
formed, unfortunately, it is unable to identify the quan- 
tum transition from the spin-singlet state (screened im- 
purity spin) to the spin-doublet state (unscreened im- 
purity spin). The slave-boson mean filed solution al- 
ways leads to a lower energy in the spin-singlet state 



than in the spin-doublet state [g6|. One must go be- 
yond the static mean-field level and employ alternatively 
the non-crossing approximation to describe this transi- 
tion (2) Very recently, Polkovnikov et al. have 
studied a similar problem. The difference and similarity 
between two works are as follows. The difference: (i) In 
their work, the Kondo exchange model where the impu- 
rity spin is coupled to the spin of conduction electrons. 
These two models should be equivalent, but the meth- 
ods to solve them are different; (ii) In addition to the 
first nearest-neighbor hopping (f), second {t') and the 
third (t") nearest-neighbor hopping were also introduced 
in their model. Instead only the first nearest-neighbor 
hopping is used in the present work; (iii) In their work, 
the modification of the superconductivity by the Kondo 
impurity was neglected, while in the present work, this 
effect has been taken into account by our self-consistent 
calculation; (iv) A vanishing density of states near zero 
bias at the sites nearest neighboring to the impurity was 
obtained in their work, while in the present work, the 
intensity of density of states on these sites was found 
to be significant. The similarity: (i) Both works have 
shown the split peaks for the local density of states di- 
rectly on the impurity site; (ii) The models in both pa- 
pers are placing the magnetic moment at the impurity 
site. This is different from the experimental situation 
for high-Tc superconductors with Zn impurities where 
the magnetic spin would most likely be induced at sites 
nearest-neighboring to the Zn impurity. Therefore, these 
two models might be more suitable to describe the effects 
of Kondo impurities such as Mn, Fe, Co substituted for 
Cu in high-Tc cuprates. 

In summary, we have considered the problem of an 
Anderson impurity coupled to the conduction electrons 
of a d-wave superconductor. Within a slave-boson mean- 
field approach, we have been able to exactly diagonal- 
ize the composite Hamiltonian with the pair breaking 
effect taken into account. A critical coupling strength, 
above which the Kondo phase occurs, is found regardless 
of whether the PH symmetry is present. In the Kondo 
effect regime, the double-peak structure appears at both 
the impurity site and its nearest neighboring sites. Our 
analysis seems to indicate that the spatial pattern of the 
differential tunneling conductance around a Zn impurity 
in BSCCO ||] cannot be fully explained with the Kondo 
impurity model. 
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